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Dynamical bistability in quantum dot structures: The role of Auger processes 
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Bistability in quantum dot structures is examined by a drift-diffusion model in combination with 
electron capture and emission processes. Our simulations provide a dynamic scenario with extremely 
long switching times of the order of months and the results are in good agreement with the exper- 
imental findings of Yusa and Sakaki [Appl. Phys. Lett. 70, 345 (1997)]. The analysis of the data 
supports the importance of Auger capture processes for quantum dots. 
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I. INTRODUCTION 

Semiconductor quantum dots (QDs) formed by self- 
organized growth processes [Q allow for the spatial con- 
finement of electrons on a scale of 10 nm. This suggests 
that these structures may serve as extremely small mem- 
ory devices. A first step in this direction is the investi- 
gation of specially grown structures where the charging 
of a layer of QDs could be detected via its influence on 
the conductivity of a two dimensional electron gas ^ |) . 
These experiments showed a pronounced bistability (up 
to room temperature 0) depending on the history of the 
bias sweep. The sensitivity to light illumination suggests 
a variety of applications in photoelectronic devices [& |) . 

A key issue towards a detailed understanding of these 
experiments is the question if the observed bistability is 
of transient nature, and the order of the time scales in- 
volved. In order to elucidate this point we performed 
a simulation of the experiment from Ref. |^ applying 
a drift-diffusion model combined with the generation- 
recombination (GR) kinetics of the QDs. Our results 
show that the experiments can be well described if Auger 
processes constitute the dominant GR-processes. There- 
fore this type of experiments may serve as an additional 
tool to shed light into the controversial issue on the na- 
ture of electron capture in quantum dots |?], ||, |^| . 



II. THE MODEL 



Here we consider the structure used in the experiment 
of Ref. || as sketched in Fig. [l| From the experimen- 
tally determined size of the QDs (height h = 5 nm and 
diameter 20 nm) we estimate a ground state binding 
energy of E h — 250 meV Q. We use a QD density 
N QD = 7.5 x 10 10 /cm 2 . 

We model the transport in the conduction band by a 
drift-diffusion approach assuming a stationary state simi- 
lar to Ref. O, which is justified if the dielectric relaxation 
time is short compared to other dynamical features. This 
provides us with the one-dimensional continuity equation 
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FIG. 1: Sketch of the sample used in Ref. 2] together with the 
resulting conduction band profile at zero bias. 
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with the characteristic function 



Xqd(z) := 



l/h (z e QD layer) 
(z <£ QD layer) 
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and the electron capture rate / which determines the 
electron density in the QD layer (in units sec _1 cm -2 ) 



i) 



Q-n&(t) = f[n(z QD ),n$ D (t)] 



(3) 



Here, e > is the elementary charge, fj, = 2 X 10 5 cm 2 /Vs 
is the electron mobility (at 77 K) and tiqq is the density 
per unit area of the electrons trapped within the QDs 
at position zqd, Ep is the Quasi- Fermi level of the elec- 
trons, which is related to the free electron density n by 



n(z) = N c (z)F 1/2 



E F (z) - E c0 (z) + e<j>{z) 



k B T 



with N r = 2 



m e k B T 
2irh 2 



3/2 
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and the Fermi integral Fi/ 2 - Here m e is the effective 
mass, E c q the intrinsic band edge of the conduction band 
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(we set E c q = for GaAs and use E c q = 0.198 meV for 
the Alg.25Gao.75As layer), T the lattice temperature, k B 
Boltzmann's constant, and (j> denotes the electric poten- 
tial. Quantization effects in the two dimensional electron 
gas are neglected for simplicity. 

The electron capture rate / is assumed to be a com- 
bination of Auger processes and single electron processes 
(like multi-phonon capture and emission) in the following 
way: 



f[n(z Q n),n™} = (Cn 2D + a){n m p^ - n 



2D„2D 
QD 



) (5) 



Here n 2B = n{zqp)h is the effective free electron den- 
sity at the QDs per unit area (which may be related 
to a wetting layer density of the same order of magni- 
tude), Pq 1 ^, = 2Aq D — Tigp is the density of unoccupied 
QD states for twofold degeneracy. C and a are the rate 
coefficients for the Auger and multi-phonon process, re- 
spectively (see also Ref. ||). For the Auger process, elec- 
tron capture requires the presence of two electrons in 
the vicinity of the QD, as well as an empty dot state. 
Thus the capture rate is proportional to (« 2D ) 2 Pqd- In 
contrast, for single electron processes the energy is trans- 
fered to different degrees of freedom (such as phonons) 
and the rate is proportional to n 2T> pqp. The respective 
electron emission rates (with negative sign) are obtained 
from the principle of detailed balance for equilibrium dis- 
tributions Jl2| . Assuming nondegeneracy in the conduc- 
tion band, we find n 2D = N c hexp[(Ep + efy/kgT] and 
"qd/Pqd — ex P[(£-F + E b + e4>)/k B T] in equilibrium. 
Thus n 2 ^/^ = N c hexp(-E b /k B T) = n 2D , and 
the prefactor n 2D in Eq. (||) ensures the vanishing of the 
rates in thermodynamic equilibrium. 

The continuity equation (|l]) has to be combined with 
the one-dimensional Poisson equation 
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with p(z) = e{Np(z) — n(z) 



dz 



-Piz) 

n QE> WXQd(z)] 



(6) 



where is the density of ionized donors resulting from 
the (5-doping and eo and e r are the absolute and relative 
permittivity. 

The boundary condition for the quasi-Fermi level is 
Ep = at the AlGaAs/GaAs interface due the con- 
tact with the two dimensional electron gas (2DEG) and 
Ep = —eV g at the GaAs/Al interface. Here the Schot- 
tky barrier height Es = 1.0 eV determines the difference 
between the conduction band edge and the Fermi level 
in the metal. Thus — e</>(zschottky) = Es — eV g provides a 
Dirichlct boundary condition for (f>. The second bound- 
ary condition of Neumann type results from the charge 
of the ionized donors at the 5-doping. All calculations 
are performed at 77 K. 
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FIG. 2: (a) Channel electron density N s versus bias V g ap- 
plied to the gate contact for a bias sweep. The duration 
of the entire sweep is one hour and an Auger coefficient 
C = 2 x 10 -12 cm 4 /s is used, (b) corresponding occupation 
probability of the QDs. 



III. RESULTS 

At first we consider only Auger processes. We use 
C = 2 x 10 -12 cm 4 /s, which was calculated in Ref. [lj| 
Fig. ^(a) shows the density N s of the 2DEG versus the 
gate voltage for a bias sweep from V g = V to 1 V and 
back to V (indicated by the arrows) with a constant 
sweep rate and a total sweep duration of 1 hour [fl4|| . 
The findings are in excellent agreement with the corre- 
sponding experimental data, see the broken line of Fig. 
1 in Ref. |^, except for an overall shift of the densities 
by ~ 2 x 10 10 /cm 2 (larger density in the experiment). 
Fig. ^(b) shows the corresponding occupation probabil- 
ity of the QDs, demonstrating that the charging state of 
the QDs discriminates between the different conducting 
states. For a comparison we have also shown the result 
for the stationary situation f[n, n 2 ^] = 0. This shows 
that both during the up- and down-sweep nonequilib- 
rium distributions persist, as the GR rates are not fast 
enough to establish thermodynamic equilibrium. 

The dynamical nature of the bistability manifests itself 
in the time dependence of the bias sweep, see Fig. ||. 
This shows that the behavior is robust against changes 
in the time scales by several orders of magnitude albeit 
the plateau density where the dots are charged is slightly 
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FIG. 3: iVs-Vg-characteristic calculated for different sweep 
durations. An Auger coefficient C = 2 x 10~ 12 cm 4 /s is used. 

changing with the sweep duration. Below this density the 
charge on the QDs is hardly changed even for the slowest 
sweep rate. This result indicates that the branches are 
stable on rather long time scales, an important feature 
for possible applications. Time scales of 100 h have been 
demonstrated in Ref. [Q at room temperature where a 
more complicated structure was used. 

The rather long relaxation times are due to the fact 
that Auger rates depend quadratically on the electron 
density n, see Eq. (||). Indeed the potential profile pro- 
vides rather low electron densities at the QDs for biases 
below « 0.7V^, leading to a very small capture rate. 

The N s (V g ) characteristic can be understood as fol- 
lows: Starting from equilibrium at V g = 0, the quantum 
dots are essentially unoccupied. This implies a linear 
potential distribution over the entire GaAs layer with a 
potential 4>(zqu) ~ (V g - E s /e)W/L at the QDs. Here L 
is the distance between the 2DEG and the Schottky con- 
tact and W is the distance between the 2DEG and the 
QDs. (We used e</>(z 2 DEG) ~ as the Fermi level Ep = 
is close to the band edge at the 2DEG.) Assuming that 
the electrons are in equilibrium with the 2DEG, we have 

For biases below V CI it this density is so small that the cor- 
responding capture rates hardly change the occupation of 
the QDs on the time scale of the bias sweep. Therefore 
the linear potential drop in the GaAs layer persists and 
the density of the 2DEG increases with V g . At V g = V CT it 
the capture rate becomes 

Then the charge captured in the quantum dots just 
screens any further increase of the bias, so that the field 
on the left side of the dots (as well as the density of the 
2DEG) remains constant. Sweeping down the bias, the 
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FIG. 4: JV s -V s -characteristic calculated with different values 
of the coefficient C with a = (full and dashed line) as well 
as C — 10 _8 cm 4 /s and a = lcm 2 /s (dash-dotted line). The 
entire sweep duration is one hour. 



electron density n{zqp) decreases and so do the capture 
rates. Therefore the charge in the QDs remains almost 
constant during the down sweep, explaining the linear 
decrease of the 2DEG electron density. 

The magnitude of the Auger coefficient is controver- 
sial. E.g., in Ref. [| a value of C = 10~ 8 cm 4 /s was 
deduced from experimental data. The corresponding 
N s (V g ) curve is given in Fig. |J (dashed curve). The result 
differs only slightly but agrees less well with the experi- 
mental data. This indicates that the findings are not very 
sensitive to the magnitude of C which just rescales the 
time scales. [This follows directly from Eq. (||): Only the 
ratio between the capture coefficients and c\V g /dt mat- 
ters.] The analysis of Ref. ^provides additionally a single 
electron capture coefficient a — Icm 2 /s in Eq. (|5|). The 
corresponding result is shown by the dash-dotted line in 
Fig. [|. The presence of the single electron term strongly 
changes the behavior. Regarding the large plateau width 
and its low electron density, this result is not compatible 
with the experiment |^]. 

In Fig. (|^) the result is given for different single elec- 
tron capture coefficients a while Auger processes are ne- 
glected (i.e. C — 0). Good agreement with the experi- 
ment is found for a = 2 x I0~ 8 cm 2 /s while values above 
a ~ I0~ 4 cm 2 /s are in strong disagreement with the ex- 
perimental findings. 

Experimentally, electron capture times of the order of 
10 — 100 ps for high carrier densities have been observed, 
see, e.g. Ref. M and references cited therein. For dot 
densities of I0 1 — 10 11 /cm 2 this implies a 3> 10 _2 cm 2 /s 
if single electron processes are dominating. From our 
study of the dynamical bistability this seems to be too 
large. On the other hand the Auger coefficient C = 2 x 
10~ 12 cm 4 /s gives capture times of 25 ps for a (wetting) 
layer density of 10 11 /cm 2 and a dot density of 10 11 /cm 2 . 
The times become even shorter for larger values of C < 
10 _8 cm 4 /s, which are also consistent with our study. 
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FIG. 5: iV s -V 9 -characteristic calculated with different values 
for the coefficient a and C = 0. The duration of the entire 
sweep is one hour. 



IV. CONCLUSION 

We have shown that the bistability observed in quan- 
tum dot structures is of dynamical nature with time 
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